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Lagrangian motions of fluid particles in a general velocity field oscillating in time are 
studied with the use of the two-timing method. Our aims are: (i) to calculate system¬ 
atically the most general and practically usable asymptotic solutions; (ii) to check the 
limits of applicability of the two-timing method by calculating the averaged motion with¬ 
out making any assumptions; (iii) to classify various drift motions and find their limits 
of applicability; (iv) to introduce a logical order into the area under consideration; (v) 
to open the gate for application of the same ideas to the studying more complex sys¬ 
tems. Our approach to study a drift is rather unusual: instead of solving the ODE for 
trajectories we consider a hyperbolic PDE for a scalar lagrangian field a{x,t), trajec¬ 
tories represent the characteristics curves for this PDE. It leads us to purely eulerian 
description of lagrangian motion, that greatly simplifies the calculations. There are two 
small scaling parameters in the problem: a ratio of two time-scales and a dimensionless 
vibrational amplitude. It leads us to the sequence of problems for distinguished limits. 
We have considered four distinguished limits which differ from each other by the scale of 
slow time. We have shown that each distinguished limit produces an infinite number of 
solutions for a{x,t). A classical drift appears in the case of a purely oscillating flow. At 
the same time, we have shown that the concept of drift motion is not sufficient for the 
description of lagrangian dynamics and some ‘diffusion’ terms do appear. Eive examples 
of different options of drifts and ‘diffusion’ are given. 


1. Introduction 

It is well-known that in oscillating flows the motion of a material (fluid) particle con¬ 
sists of two parts: oscillating and non-oscillating. The latter is known as a drift, see 


Stokes (1847)1 Maxwell (1870)[ Lamb (1932)] Longuet-Higgins (1953)[ Darwin (1953)[ Hunt (1964) 


Batchelor (1967)[Lighthill (1956)[ Andrews fc McIntyre (1978)[ Craik (1982), Grimshaw (1984)[ 


Craik (1985) Benjamin (1986)[ Eames, Belcher, and Hunt (1994), Eames fc McIntyre (1999)| ~ 

Buhler (2009)[ and many others. The great significance of a drift due to its key role in 
the vortex dynamics of oscillating flows and Langmuir circulations was recently analyzed 


from a new perspective by Vladimirov (2012) Vladimirov, Proctor and Hughes (2015) 


Here we study the advection of a lagrangian scalar field by a given oscillating velocity. We 
employ the two-timing approach combined with eulerian time-averaging operation in the 
form introduced by Vladimirov (2005) Yudovich (2006) [ Vladimirov (2008) The two- 
timing method has been used by many authors, see Nayfeh (1973) Sanders fc Verhulst (198^ 
however our analysis goes well beyond the usual calculations of main approximations in 
various special cases. Our analytic calculations are straightforward by their nature, but 
they do include a large number of integration by parts and algebraic transformations; 
the performing of such calculations in a general formulation represents a ‘champion-type’ 
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result by itself. These calculations are too bulky to be presented in JFM, they are de¬ 
scribed in detail in the arXiv papers by Vladimirov (2010)] Vladimirov (2011) quoted 
below as I and II, therefore here we present only some results of these calculations. The 
main purpose of this paper is to introduce a new systematic and general viewpoint on the 
hierarchy of distinguished limits and drifts. This new viewpoint helps to organize and 
unify a number of results, including that of I and II. 


2. Formulation of Problem 

A fluid flow is given by its velocity field v*(x*,t*), where x* = {x\,X 2 ,x’^) and t* are 
cartesian coordinates and time, asterisks stand for dimensional variables. We suppose 
that this field is sufficiently smooth, but we do not suppose that it satisfies any equations 
of motion. The dimensional advection equation for a scalar lagrangian field a{x*,t*) is 

Of-t (u* • V*)a = 0, af^=da/dt* (2.1) 

This equation describes the motions of a lagrangian marker in either an incompressible 
or compressible fluid or the advection of a passive scalar admixture with concentration 
a{x*,t*) in an incompressible fluid. In the case of a compressible fluid it also describes 
the advection of a passive scalar admixture, where a represents a ratio of concentration of 
the admixture to the density of a fluid. The hyperbolic equation has characteristics 
curves (trajectories) x* = x*{t*) described by an ODE 

dx*/dt*=v\ x*\t=Q = X* (2.2) 

where x* and X* are eulerian and lagrangian coordinates. The classical description of 
drift motion follows after the integration of (12.2p . however for higher approximations it 
requires very bulky operations with lagrangian displacements. Therefore we solve the 
equation (ED that allows to use the eulerian average operation and to simplify calcula¬ 
tions. The oscillating field v* possesses the characteristic scales of velocity U, length L 
and frequency w*. These three parameters give a Strouhal number St = Loj*/U, hence 
the dimensionless variables and parameters (written without asterisks) are not unique; 
we use the following set 

X* = Lx, t* = {L/U)t, uj* = {U/L)bj, V* = Uoj^u = {St)^Uu; ft = const < X2-3) 

We accept that u has the ‘two-timing’ functional form 

u{x,s,t), t = wt, s = t/uj°‘; M ^ 0(1), a = const > — 1 (2.4) 

where r-dependence is always 27r-periodic, while s-dependence is arbitrary. Two indefinite 
constants a and ft will be defined in distinguished limits. The restriction a > — 1 makes 
the variable s ‘slow’ in comparison with r, while ft < \ gives a small vibrational spatial 
amplitude; uj = St is considered as a large parameter. 

In dimensionless variables (after the use of chain rule) (12.11) takes the form 

wOt H- -Os + [u ■S/)a = Q, d/dt = ujd/dr +-^d/ds (2.5) 

where the subscripts r and s stand for partial derivatives; s and r represent two mutu¬ 
ally dependent time-variables, which are called slow time and fast time. The functional 
structure of u (j2.4l) underpins the key suggestion that the solution of (12.11) has the same 
structure 


a = a{x, s, t) 


( 2 . 6 ) 
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Equation (I2.5|l can be rewritten in the form containing two independent small parameters 
Or + £20^ + ElM • Va = 0; El = £2 = (2.7) 


Hence, we operate in the plane (£i,£i), where we study asymptotic limit (£i,£ 2 ) —>■ 
(0,0). Different asymptotic paths (£i,£ 2 ) —t (0,0) may produce different solutions, such 
paths are called distinguished limits. In (12.71) . following the methodology of two-timing 
approach, we (temporary) consider r and s as independent variables. 

In order to make further analytic progress we require some specific notations. In this 
paper we assume that any dimensionless function f{x, s, r) has the following properties: 
(i) / ~ 0(1) and all its required x-, s-, and r-derivatives are also 0(1); (ii) is 27r-periodic 
in T, i.e. f{x, s, t) = f{x, s,t + 27r); (iii) has an average given by 


(/)^ 


f*To+27r 


27r, 


f{x,s,T)dT = f{x,s) 


V To = const; 


(iii) can be split into averaged and purely oscillating parts, f{x, s, r) = f{x, s) + f{x, s, r) 
where the tilde-functions (or purely oscillating functions) are such that (/) = 0 and the 
bar-functions f{x, s) are r-independent; (iv) we introduce a special notation for the 
tilde-integration, which keeps the result of integration in the tilde-class: 

r = /(®,s,T')dT '-(^J f{x,s,T')dT''jdfi. (2.8) 


3. Distinguished Limits and Related Solutions 

It can be shown (see I,II) that there are series of distinguished limits for the equation 
(lO) . which represent the one-parametric distinguished path £i = £, £2 = £" with an 
integer n > 0: 

Or-\-e’^Qs-\-eu ■ \7a = 0 (3.1) 

Let us consider first four cases of Distinguished Limits DL(n) with n = 1,2, 3,4: 

DL(1) s = t] u = u(x, s) -\- u{x, s, t); 

DL(2) s = et] u = u{x, s, r); 

DL(3) s = £^t; u = u{x,s,t) which satisfies condition Vq = 0; 

DL(4) s = £^t; u = u{x,s,t) which satisfies conditions Vo = 0, Vi = 0. 

The used notations are 

Vo = i[S,i], Vi = i([[S,|],|]), | = (3.2) 

where the commutator [/, g] = {g ■ V)/ — (/ • V)g for any vector-functions / and g. In 
all cases we are looking for the solution in the form of regular series 

00 

a{x,t,T) ='^e''ak{x,t,T), fe = 0,1,2,... (3.3) 

fc =0 

The case DL(1) with u = 0(1) ^ 0 naturally corresponds to the advection speed of 
order one, hence s = t and the averaged equation of zero approximation is 

(ds + M • V)ao = 0. 

Here we do not consider here this case in detail, one can find it in 1. 

The most instructive from our point of view is the case of DL(2). Here we have a longer 
slow time scale s = et, since here u = 0 and, naturally, the speed of advection is lower. 
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The substitution of (I3.3p into (131]) and a number of transformations lead to the averaged 


equations of successive approximations, which are derived in I, II 

{ds + Fo • V) ao = 0, ds = d/ds (3.4) 

{ds + Fo • V) ai + {Vi ■ V)ao = 0 (3.5) 

(a« +Vo-V)a2 + (Vi • V)ai + (^2 • V)ao = ^ > (3-6) 

^2 = ^([[^0,?],!]) + \{[Vo,vl]) + + i(idivM +MdivO, (3.7) 

S' = |,-[Yo,e (3.8) 

‘^Xik = 'Vu'k^i) = j (3-9) 

%J.fc = {ds + Fo • V) 7,, - ^l.ra - ^7km (3-10) 

UXm 


where the operator is such that £'Ygfik — 9 represents the condition for tensorial 
field /jj,(a;,t) to be ‘frozen’ into Vo{x,t) (also known as the Lie derivative). Three equa¬ 
tions (isa-drel) can be written as a single advection-pseudodiffusion equation (with an 
error 0(e^)) 

V = V'"’= Vo + eVi+e2y2, (3.12) 

Kzk=xfk=£%k (3-13) 

a = =00-1- eoi -I- £^02 (3.14) 

_ _| 2 ] 

Eqn. (13.1111 shows that the averaged motion represents a drift with velocity V = V + 
0(£^) and pseudodiffusion with matrix coefficients 'Kik = s'^Xik + 0(e^^)- We have intro¬ 
duced the term pseudodiffusion on the ground of the following: (i) the evolution of a 
is described by an advection-diffusion-type equation (13.111) : (ii) the diffusion type term 
represents a known source-type term in the second approximation; and (iii) the equation 
(|3.1ip is valid only for regular asymptotic expansions (I3.12ll - (l3.14p . 

For the case DL(3) we impose a restriction Vq = 0 and derive equations (see I): 


(a.-kVi-V)ao = 0 (3.15) 

(a. + • V) ai + (F 2 • V)ao = ^ , (3-16) 

V 2 ^ i([[Vo,?],!]) + i([Vo, + ia.(idiv0, (3.17) 

2x.fc = dsiMk) (3.18) 

For the case DL(4) we impose two restrictions Vo = 0 and Vi = 0 and derive the 
equation (see I): 

(a.+V..V)S. = |-(&g), (3.19) 


with the same V 2 and Xik as in DL3. The comparison between the averaged equations 
for DL(1)-DL(4) shows that the same coefficients from higher approximations in DL(N) 
appear in DL(N-l) in the equations of the previous order in e. The higher approxima- 
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tions DL(5), etc.) can be derived similarly, however the calculations become extremely 
cumbersome. 

One should recognize, that in all presented distinguished limit solutions the meaning 
of key parameter e has not been fully defined yet. To make a progress in this direction, 
let us again consider DL(2) for (13.11) 

Or + + eu ■ Va = 0 (3.20) 

In order to make this equation coinciding with (12.51) we must choose a and /3 as: 

/?=(l-a)/2, e = (3.21) 

Now, for any solution of (I3.20L (which represents a distinguished limit solution) we 
obtain an infinite number of solutions, one solution for any real number a > —1. Those 
solutions correspond to different slow time variables and different magnitudes of given 
velocity. For example: 

• If M is given as function of variables t = ujt and s = t, then we must take a = 0, 
/3 = 1/2, hence a dimensionless velocity {St)^u in (12.31) is of order 0{^/e), and e = \fiw. 

• The most frequently considered case, when a velocity is of order one, corresponds to 
a = 1, /? = 0, and e = 1/w. 

• Another interesting possibility corresponds to the case a = /3 = 1/3. Here, s = tf 

e = and velocity {St)^u = Although such a scaling may look exotic, it is 

required if a particular slow time-scale s = tj is prescribed in u. 

• Our results for DL(2) are particularly striking for the case of u{x,t) not depending 
on s, which represents a flow with purely periodic velocity oscillations. As one can see 
from the above consideration, the scale of slow time in this case is uniquely determined 
by given velocity lo^u. For every value of /3 we have the related slow time s = tfu}'^. It 
must be accepted, that solutions with different /3 are physically different. 

• The general tendency given by (13.211) is physically natural: in order to shorten the 
slow time-scale (decreasing a), one needs to increase the amplitude of a given velocity 
(increasing j3). 

• In general, transformations similar to (I3.20I1 . (I3.21|) produce an infinite number of 
different solutions for each case DL(l-4). 


4. Examples 

Example 1. The superposition of two modulated oscillatory fields of the same frequency. 

u{x, s, t) = p(x, s) sinr -|- q(x, s) cost (4.1) 

where p and q are arbitrary bar-functions. The straightforward calculations yield [m, = 

\p,q]. The drift velocities (13.2|) . (13.7p are 

Vo = ^[p,ql Ei = 0, (4.2) 

^2 = i {[P,p] + [Q,q]) - 2 ([Ps>p] + [Qs^Q]) + (4-3) 

+ ^ ^pdiv P' + qdiv q' + P'divp + Q'divq'j , 

P=\Vo,pi Q=\Vo,ql P' = Pg-P, Q' = q,-Q, 

(Mk) = -^{PiPk+Wk) 


(4.4) 
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Pseudo-diffusion matrix Kik follows after the substitution of (14.41) into (13.91) . The velocity 
u (14.1[) is general enough to produce any given function Vo{x,t). To obtain p{x,t) and 
q{x,t) one has to solve the equation 

[p,q] = (q • V)p- {p =2Vo{x,s) (4.5) 


which represents an undetermined bi-linear PDE-problem for two unknown functions. 
Example 2. Stokes drift. 


The dimensionless plane velocity field is (see Stokes (1847) Lamb (1932) Debnath (1994)) 
u = Ae^y 


cos{kx — t) 
sin(fca; — r) 


{x,y) = {xi,X 2 ) 


(4.6) 


where one can choose A = 1 and fc = 1; however, let us keep both A and k in the formulae 
in order to trace the physical meaning. The fields p(x, y), q{x, y) (14.ip are 


— A kv f smkx 
p = Ae^y ( , 

— cos kx 

The calculations of (14.21) yield 

Vo = kA^e^’^y 


Q = 


. ky ( coskx 
Y sin/ex 


Vi = 0 


(4.7) 


(4.8) 


which represents the classical Stokes drift and a zero first correction to it. For brevity, 
the explicit formula for V 2 is not given here. Further calculations show that 

Xik = -X , with X = 

One can see that the eigenvalues Xi = ~X X 2 — X correspond to a strongly 
anisotropic ‘diffusion’. The averaged equation (13.111) (with an error 0{e^)) can be written 
as 


at + {yo-\-e‘^V2)aj: = e^{Xyax + xaxy) (4.9) 

(X = do “1“ £ai “h 0,2 


where V o and V 2 are the x-components of corresponding velocities (their y-components 
vanish). This equation has an exact solution a = a{y) where a[y) id an arbitrary function, 
which is not effected by pseudodiffusion. 

Example 3. A spherical ‘acoustic’ wave. 

A velocity potential for an outgoing spherical wave is 

~ A 

(j) = — sin(fcr — r) (4.10) 

r 

where A, fc, and r are amplitude, wavenumber, and radius in a spherical coordinate 
system. The velocity is purely radial and has a form dUD 


u = p sin T + q cos r, 
p = A \ — cos kr — sm kr 

\ J.2 j. 


q = A \ —- sin kr -\— cos kr 


(4.11) 

(4.12) 


whCTe u^p, and q are radial components of corresponding vector-fields. The fields ^ and 
[u, are also purely radial; the radial component for the commutator is 


^Ur — U^r = A'^k^/r'^ 


(4.13) 
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where ^ is radial component of ^ and subscript r stands for the radial derivative. The 
drift ()4.2I) is purely radial with 


Vo = 




^1 = 0 , ^2 = 


16 r^ 


3r - 


(4.14) 


It is interesting that V o formally coincides with the velocity, caused by a point source in 
an incompressible fluid, and for small r the value of 1^2 dominates over I/qi so the series 
is likely to be diverging. Further calculations yield 

= + X = 4 l 4 fcV 4 r 2>0 ( 4 . 15 ) 

where y stands for the only nonzero rr-component of Xife. One can see that in this case 
pseudodiffusion appears as ordinary diffusion. 

Example 4- Vi-drift. 

It is interesting to see such flows for which Vo = 0 but Vi 7 ^ 0. Let a velocity field be 
a superposition of two standing waves of frequencies lo and 2 a;: 

u(x, t, t) = p{x, t) sinr + q{x, t) sin 2t (4-16) 

[u,$,] = i [p, q] (2 cos r sin 2r — cos 2t sin r) (4.17) 

Hence (13.21) yields 

^0 = i([M,l]) = 0, Vi = = g[[p,q],p] (4.18) 

These expressions produce infinitely many examples of flows with Vi-drift. 

Example 5 Chaotic dynamics for V 0 . 

Let a solenoidal/incompressible velocity (14.11) be 


P = 


cosy \ 

j 

f a sin 2 

0 , 


b sin X -I- a cos z 

siny / 

1 

[ b cos X 


(4.19) 


where (x, y, z) are Cartesian coordinates; a, b are constants. Either of this fields, taken 
separately, produces simple integrable dynamics of particles. Calculations yield 

( —asinysinx — 26sinycosz \ 

b sin 2 ; sin y — a cos X cos y , (4.20) 

b cos 2 : cos y + 2 a sin x cos y J 


Computations of solutions for dxo/ds = Vq which is following from (12.21) for DL(2) (see 
I) were performed by Prof. A.B.Morgulis (private communications). He has shown that 
this steady averaged flow exhibits chaotic dynamics of particles. In particular, positive 
Lyapunov exponents have been observed. Hence, the drift created by a simple oscillatory 
field can produce complex averaged lagrangian dynamics. 


5. Discussion 

1. Our study is aimed to create a general viewpoint on distinguished limits and drifts 
based on the two-timing method. We do not use any additional suggestion and assump¬ 
tions, hence any physical failure of our results can be related to the insufficiency of two 
time-scales only. Indeed, the presence of scaling parameters, such as St, ei, and £ 2 , allows 
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to introduce an infinite number of additional time-scales. From this perspective the prob¬ 
lem we study can be seen as a test for sufficiency of the two-timing method. In particular, 
some secular (in s) terms could appear in solutions of (13.61) and (13.161) . If it is recognised 
as unacceptable, then further time-scales should be introduced, which require the sys¬ 
tematic development of three-timing etc. methods. Alternatively, one can suggest that 
the two-timing method fails at the orders of approximations, where the pseudo-diffusion 
appears. 

2. One can rewrite the approximate solution (13.141) back to variable t (along with its 
tilde-part, see I and II) and substitute it back into the exact original equation (12.1|) : then 
a small residual (a nonzero right-hand-side) appears. The method used allows to produce 
an approximate solution with a residual as small as needed. However, the next logical 
step is more challenging: one has to prove that a solution with a small residual is close to 


the exact one. Such proofs had been performed by Simonenko (1972) Levenshtam (1996) 


for the case of vibrational convection. Similar justifications for other oscillatory flows are 
not available yet. 

3. All results of this paper have been obtained for the class of r-periodic functions, 
which is self-consistent. One can consider more general classes of quasi-periodic, non¬ 
periodic, or chaotic solutions. The discussion on this topic is given in I, II. It is worth to 
understand the properties of periodic oscillations, in order not to link these properties 
exclusively to more general solutions. From other side, the simplicity of calculations for 
the T-periodic solutions allows to obtain very advanced results that can be searched and 
generalized for more general solutions. 

4. The appearance of pseudo-diffusion in the DL(I-4) is especially interesting to study. 
Its physical mechanism is discussed in detail in I, II. It is worth to mention that the 
averaged pseudo-diffusion coefficient for different flows can be positive, negative, or it 
can change sign in space and time. 

5. Example 5 demonstrates that a drift can produce chaotic averaged dynamics of par¬ 
ticles. This result leads to numerous new questions and opportunities such as: (i) what 
is the relationship between chaotic motions for an original dynamical system and for the 
averaged one? (ii) how can a chaotic drift and pseudo-diffusion complement each other? 
(iii) how a chaotic drift can be used in the theory of mixing? (iv) since the averaged dy¬ 
namics is chaotic, then the related results by Arnold (1964)[ Aref (19M)| Ottino (1989) 


Samelson & Wiggins (2006) Chierchia fc Gallavotti (1994)|can be efficiently used. 


6 . It is worth to calculate the characteristics/tragectories directly by solving (12.21) with 
the use of the same two-timing method and to compare the related solutions for drifts 
with the results presented above. Such calculations are partially presented in I, II. 

7. Similar results for a passive vectorial admixture are presented in I, II. They are 


closely linked to the problem of kinematic MifU-dynamo, see Moffatt (1978) It is phys¬ 
ically apparent, that for the majority of shear drift velocities V{x,t) the stretching of 
material elements produces the linear growth of a magnetic field |ft,| ~ t. At the same 
time, for the averaged flows with exponential stretching of ‘averaged’ material lines these 
examples will inevitably show the exponential growth. 


The author is grateful to Profs. A.D.D.Craik, R.Grimshaw, K.I.Ilin, M.E.McIntyre, 
H.K.Moffatt, A.B.Morgulis, T.J.Pedley, M.R.E.Proctor, and D.W.Hughes for helpful dis¬ 
cussions. 
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